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Abstract 

We study the production and equilibration of the gluon-minijet plasma ex- 
pected to be formed in the central region of ultrarelativistic heavy-ion col- 
lisions at the BNL-RHIC and the CERN-LHC by solving a self-consistent 
relativistic transport equation. We compute the minijet production within 
perturbative QCD. Subsequent collisions among the semi-hard partons are 
treated by considering the elastic gg — > gg processes with screening of the 
long wavelength modes taken into account. We determine the time Teq where 
close to ideal hydrodynamic flow sets in, and find rather similar numbers 
for central heavy-ion collisions at BNL-RHIC and CERN-LHC energies. The 
number densities, energy densities and temperatures of the minijet plasma are 
found to be different at RHIC and LHC, e.g. T{Teq) ~ 220 MeV and 380 MeV, 
respectively. 

PACS: 12.38.Mh, 25.75.-q, 13.87.-a, 24.85.-hp 
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I. INTRODUCTION 



In the next few years, the BNL-RHIC (Au-Au colhsions at ^/s=200 GeV per incident 
nucleon pair) and the CERN-LHC (Pb-Pb colhsions at 1/5=5. 5A TeV) accelerators will 
provide the opportunity to study a new phase of matter, namely the so-called Quark-Gluon 
Plasma (QGP) fl]]. It is very important and interesting to study whether the QGP actually 
does thermalize in those reactions, and if so, what is the actual energy density, number 
density and temperature at which it thermalized. For this purpose, and also for the calcu- 
lation of all the signatures, it is necessary to study the space-time evolution of partons just 
after the nuclear collision. For example, the equilibration time is crucial for a quantitative 
understanding of J/ip suppression [0, and it is a challenging task to determine this quan- 
tity accurately. Similarly, understanding the equilibration time is important for all other 
proposed signatures for the QGP, for example for dilepton emission and strangeness pro- 
duction Once equilibrium is reached, the further space time evolution of partons can be 
described by the well known equations of hydrodynamics. 

The evolution of the QGP towards (local) equilibrium can be studied by solving transport 
equations for quarks and gluons with all the dynamical effects taken into account. Obviously, 
the first problem one always encounters is the correct computation of the initial conditions 
needed to solve the transport equation. This is because one can not calculate the parton 
production in all range of momentum from perturbative QCD (pQCD). There are also 
coherence effects that play an important role in the early stage of the nuclear collision 
at very high energy. For small x and large nuclei, the QCD based calculation performed 
in 1^ predicts the existence of a coherent field in a certain kinematical range. That field may 
play an important role in the equilibration of the plasma. In the present paper, however, 
we restrict our calculation to the initial incoherent parton production, which is computed 
within the framework of pQCD. We study the subsequent evolution of that minijet "plasma" 
by solving a relativistic transport equation, thus taking into account collisions between the 
produced partons. In the future, we intend to generalize our approach to include both 
coherent field and incoherent partons in the transport equation. 

The paper is organized as follows. In section II we briefly review minijet production 
in high-energy nuclear collisions within pQCD. In section III we discuss the in-medium 
screening of long wavelength gluons which is relevant to our study. We present the relativistic 
transport equation in section IV, briefly describing the numerical strategy for its solution in 
section V. We discuss our main results in section VI and conclude in section VII. Throughout 
the manuscript we employ natural units, h = c = k = 1. 



II. MINIJET PRODUCTION IN NUCLEAR COLLISIONS AT RHIC AND LHC 

In this section we review the computation of the single- inclusive semi-hard cross section 
in lowest order pQCD, cf. also PJ^l- "^^^ 2 — »• 2 minijet cross section per nucleon in AA 
collision is given by 

ajet = j dptdyidy2^^^^xi fi/A{xi,Pt) ^2 fj/A{x2,Pt) aij^ki{s,i,u). (1) 

^ ijkl 
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Here xi and X2 are the light-cone momentum fractions carried by the partons i and j from 
the projectile and the target, respectively. fj/A are the distribution functions of the parton 
species j within a nucleon bound in a nucleus of mass number A. yi and 1/2 denote the 
rapidities of the scattered partons. The symbols with carets refer to the parton-parton cm. 
system, dij^ki is the elementary pQCD parton cross section. 

s = X1X2S = Apl cosh^ \ ' (2) 



gives the total c.m.-energy of the parton-parton scattering. The rapidities yi, y2 and the 
momentum fractions Xi, X2 are related by, 

xi = pt {e^' + 6^')/^/^, X2 = Pt {e-y' + e-y')/V^. (3) 



The limits of integrations of rapidities yi and y2 are given by \yi\ < \\i{y/s / 2pt + \J s j Apl — 1) 
and — ln(yi/]9( — e"^^) < y2 < ^^{y/s/pt — e^^), respectively. We multiply the above minijet 
cross sections by the phenomenological factor K = 2 to account for higher-order contribu- 
tions. 

The minijet cross section, Eq. (|I]), can be related to the total number of produced partons 

via 

N = T(0) f dptdyidy2^^^xi fi/A{xi,Pt) X2 fj/A{x2,Pt) aij^ki{s,i,u), (4) 

^ ijkl 

where T(0) = dA^ /87iR\ is the nuclear geometrical factor for head-on AA collisions (for a 
nucleus with a sharp surface). Ra = 1.1 A^^^ fm is the nuclear radius. Similarly, the total 
transverse energy E of minijets is given by 



dptPtdyidy2—::^^Xi fi/A{xi,pl) X2 fj/A{x2,Pt) aij^kiis,i,u). (5) 

^ ijkl 



We employ the "EKS98" set of parton distribution functions in a bound nucleon from ref . , 
based on the GRV98 set of parton distributions for a free nucleon One has to choose the 
scale in the transverse momentum below which the incoherent parton picture is not valid 
and coherence effects have to be taken into account. These scales provide the initial cutoff 
for the (semi-)hard scatterings. The values for the cutoff are found to be ~ 1 at RHIC 
and ~ 2 GeV at LHC from the McLerran-Venugopalan model using their initial conditions 
For the quantitative estimates presented below we actually choose the cutoff po to be 
1.13 GeV and 2.13 GeV at RHIC and LHC, respectively. These values are obtained by an 



independent method by Eskola et. al, [10 



To solve the transport equation one has to specify the initial distribution function of the 
particles in phase-space, while accounting for the correlation between momentum-space ra- 
pidity and space-time rapidity. We choose our initial distribution function to be a Boltzmann 
distribution in the local rest frame, 

/(to = 1/po, i,Pt) = exp {-p^u^'/Tjet) = exp {-pt cosh^ /Tjet) ■ (6) 

m'^ = (cosh?7, 0, 0, sinh?7) is the four- velocity of the local rest-frame of the medium, with 
1] = Artanh (z/t) denoting the space-time rapidity, and C, = rj — y. The parameter Tjet can 
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be determined from the average energy per particle of the initially formed minijets which is 
obtained from Eqs. and @. We mention here that even though the initial distribution 
is chosen to be a Boltzmann distribution, the subsequent evolution will be non-ideal, see 
below. 

We shall use these initial conditions to study the evolution of the parton "plasma" 
at RHIC and LHC by solving a self-consistent relativistic transport equation (see section 
IV). As most of the produced minijets are gluons, we will simplify our considerations by 
considering gluons only. In Table-I we have listed the initial conditions of the gluon minijets. 



III. SCREENING IN NON-EQUILIBRIUM 

In this section we describe the screening of long wavelength electric fields in the parton 
medium. It will play an important role in defining the finite collision term of our transport 
equation. We shall employ the static limit (screening of infinitely long wavelength fields), 
which gives the Debye screening mass, and assume that this simplified treatment is at least 
qualitatively correct. Our concern here is to incorporate the density dependence of the 
medium- induced cutoff (which thus also depends on collision energy), as well as to treat 
the coupled evolution (cutoff and medium) self-consistently. At present we can not include 
magnetic (dynamic) screening into the non-equilibrium study because we do not find an 
expression for the magnetic screening mass in terms of the non-equilibrium distribution 
function in the literature (see section VII for a detail discussion). We also mention here 
that the momentum dependent screening mass is usually studied in thermal equilibrium. 
However, accounting for the momentum dependence of the screening mass in non-equilibrium 
is technically very difficult and is beyond the scope of this present paper. 

The electric screening mass is given by the infrared limit of the real part of the gluon 
self-energy U^^, calculated in the given background that is described by the distribution 
function / (not to be confused with the parton distribution function fj/j^ introduced in 



section II). In ref. ||Tl| the following expression has been derived (in Coulomb gauge) for a 



medium of gluonic excitations: 

^2 ^ _3«. j.^ r ^3 i V /(p). (7) 

vr^ \q\^oJ q ■ P 

In the above equation g is the momentum of the test particle, f{p) is the non-equilibrium 
distribution function of the gluons and is the strong coupling constant. We will consider 
the transverse screening mass in the following, which will be introduced below as a cutoff 
in parton-parton anti-collinear elastic scattering to obtain a finite transport cross-section. 
Performing an integration by parts we obtain the transverse screening mass 

Sag f d'^p 



r-^m. (B) 



TT^ J p 



Following Bjorken's hypothesis [Q, we express all quantities in terms of longitudinal-boost 
invariant parameters r = \/t'^ — z'^, ^ and pt- We assume that the above equation is also 
valid for a space-time dependent distribution function f{x,p) and hence use 



2/ N 6a, (r) 
(r) = — — 



JdptptJdUXr,pt,0 ■ (9) 
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Improving earlier approaches |T3[ , we do not assume factorization of the distribution function 
in the form f{T,pt,^) = g{i)h{T,pt). Also, in our case the screening mass enters the collision 
kernel of the transport equation and thus determines the rate of equilibration. 

In the above equation the QCD running coupling constant becomes time dependent, via 

a,(r)=a,((p?(r))). (10) 

The average transverse momentum squared of the excitations of the medium is defined as 

Here dT = dV(2vr)V 

= dpjPid^/(27r) is the invariant momentum-space measure. 
Throughout the manuscript we use the GRV98 calculation of as{{pf{T))) ^ with {pjij)) 



calculated from Eq. (ITT]). 



IV. SOLUTION OF THE TRANSPORT EQUATION WITH SCREENING 

In the absence of any coherent color field the space-time evolution of the produced partons 
at RHIC and LHC can be studied by solving the Boltzmann transport equation 

p^d,f{x,p) = C{x,p), (12) 

where f{x,p) is the distribution function of gluons and C{x,p) is the collision term. To 
solve the above transport equation with the initial value of /o given by Eq. (^, we employ 
the relaxation time approximation for the collision term [lT4| , |T5i : 

C(r, pt) = -p^u, [/(r, i^pt)- rir, ^, p,)] /t,{t) . (13) 

f'^'^ is the Bose-Einstein equilibrium distribution function, and Tc(r) is the time dependent 
relaxation time of the plasma. With this collision term the formal solution of the transport 
equation becomes 

/(.,J,p,)^£d.'exp(/;'-|l)^!^+ /K5,.p,)exp(-/;-|l), (14) 

where ^' is the solution of 

sinh^' = — sinh^ , (15) 

and sinh^o = :^sinh^, with Tq = I/pq. Despite the fact that the transport cross section 
(governing kinetic equilibration) has been derived from the gg — >■ gg elastic scattering cross 
section, the number-current of the gluons is not conserved. Rather, for massless particles 
the number density in equilibrium is proportional to the entropy density. However, dissi- 
pation during the pre-equilibrium stage will produce additional entropy, and accordingly 
the comoving number density is expected to decrease less fast than for a conserved current 
= pu^, where (9^j(pM^) = dp/dr + p/r = 0, see below. 
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We write the relaxation time for collisions, Tc(r), as 

1 



at{T)n{T) 



where 



n{T) = 9g J dTpf,u''f{T,pt,0 
is the number density of the gluon-minijet plasma, and 

(^tir) = / sin^ Ocm. 



(16) 



(17) 



(18) 



denotes the time dependent transport cross-section for the collision processes |T^JT^. We 
assume that anti-collinear small-angle scattering gives the dominant contribution to the 
transport cross-section |T^, such that sin^ 6'c.m. = ^tii/s^. We mention again that, in our 
study, all the quantities such as m^, at, n, Tc are time dependent and have been obtained 
from the non-equilibrium distribution function of the gluon minijets. 

We shall consider the leading order elastic scattering processes gg gg. The differential 
cross section for this process is given by 



da 
dt 



~2P 



ut 



us 



St 



u 



(19) 



In the limit of small-angle scattering (of identical particles) it simplifies to 

dcr Qvrag 1 
dt ~ 



(20) 



and the transport cross-section cXf diverges logarithmically due to exchange of long- 
wavelength gluons. However, as discussed in section III, long-wavelength fields will be 
screened by the dense medium. For our studies we therefore employ the medium-modified 



elastic cross-section [17 



da 



(21) 



Using Eq. (|21|) in Eq. ([T8|) we obtain the medium modified finite (time dependent) transport 
cross-section 



9 47r«2(^) /^2(^) 



1 



(s(r)+2m2(r)) log 





, / sir) 






-f- 1^ - 2s(r) 









2 s2(r) V Kr) 

To simplify the considerations we have replaced s by its average value 

sir) = A{E{r)f 



(22) 



(23) 



in the above expression. {E{t)) is the time dependent average energy per particle given by 
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(24) 



where 

e(r) = (7G/dr(pX)V(r,P*,0 (25) 

is the energy density of the minijet "plasma" and n(r) is the local number density as defined 
in Eq. (^). 

To obtain the collision term (|13|) in each time step we also have to determine the equi- 
librium distribution towards which the evolution is supposed to converge. In other words, 
we have to determine the "equivalent" plasma temperature from (pS]). It can be obtained 
from the condition that the first moment of the actual distribution function / equals that 
of the equilibrium distribution function Z^"^, i.e. 

This should be a reasonable approximation as long as the system is not very close to the 
hadronization phase transition. 



V. NUMERICAL SOLUTION 

The expression for the distribution function f{T,C,,Pt), Eq- (0)? involves T(r) and Tc{t) 
which are again defined through the distribution function /(r, ^,pt)- We solve these coupled 
set of equations self-consistently. At any time r we start with the old trial values To, no, 
«sO) "^to and So from which we get fo via Eq. (|14D. This fo is used in Eq. ( p5D to calculate 
e(r) which gives a new temperature T/v through Eq. (pG]). This new temperature, but the 



old values of no, ctso^ "^to ^-nd so are again used in Eq. ([T^) to get a new fi. This /i is 
used in Eq. (|r7p to calculate a new n^r which also gives a new value via Eq. (^). These 
new values T^r, U]^ and sn and the old values aso and mto are again used in Eq. (|I^ to get 
a new /2. Using this /2 in Eq. (|TTp we obtain asN from Eq. (^). These new values T/v, njsi, 
Sn, otsN and old value mto are again used in Eq. (|^) to obtain a new /s. Using this /s in 
Eq. (H) we obtain a new mt^. Thus, starting with the old set of values To, no, otsOi ^to and 
So we obtain a new set of values Tjy, n^, cXsn^ t^in and s^- This process is iterated until 
convergence is attained to the required accuracy. This gives us the self-consistent values of 
e(r), T{t), nij), as(r), mtir) and s(r) at any time r. 



VI. RESULTS AND DISCUSSIONS 

The purpose of this paper is to study various bulk properties of the gluon minijet plasma. 
In particular, we discuss the time evolution of the number density, the energy density, 
the temperature, as well as the collision-relaxation time and the transport cross section. 
Calculations of signatures from this equilibrating minijet plasma will be presented elsewhere. 

We compare the time evolution of the energy density with that obtained in the free 
streaming limit (no collisions and initial rapidity distribution dN/dy oc 6{y — rj)) and in the 
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equilibrium limit (isotropic momentum distribution in the comoving frame, at all times), i.e. 
the hydrodynamical evolution. For purely longitudinal expansion the energy density in the 
free streaming limit behaves as 

(27) 

where is the energy density of the gluon-minijet plasma at the time Tfg where free streaming 
sets in. 

In the equilibrium limit we have 

n{T)=nJ—] , (28) 




eW = eo(^^j , (29) 

T(r)=To(^-^) \ (30) 

In the above equations eo, and Tq are the energy density, number density and temperature 
of the gluon plasma at r = Tgg, where equilibrium is reached. The latter two equations 
actually depend on the equation of state of the minijet plasma; we have assumed an ideal 
gas of gluons. 

The evolution of the energy densities at RHIC and LHC are shown in Fig. 1 and Fig. 2 
respectively. The solid lines depict the result from our self-consistent transport calculations. 
The dashed lines correspond to the free streaming evolution (see Eq. (p7|)) with the minijet 
initial conditions at r/s = tq = 1/po- The dot dashed lines correspond to the 101 hydro- 
dynamical evolution of the energy densities (see Eq. ([29|)) with the same initial conditions 
of the minijet plasma at Teq = tq. The latter case represents the evolution of the energy 
densities assuming equilibration at r = tq. It can be seen that our results lie between the 
free streaming and hydrodynamic limits. To see when equilibration sets in we fit our results 
with e(r) oc r~", at different times. We present our fitted values of a in Table-II and Table- 
Ill for RHIC and LHC respectively. It can be observed that as time progresses the scaling 
exponent a increases. This is due to the fact that at later times the collision rate among 
the plasma constituents increases, driving the system towards equilibrium. As a cross-check 
we also examine the behavior of the temperatures and number densities. 

The time evolutions of the "temperatures" of the minijet plasma at RHIC and LHC are 
depicted in Fig. 3. It can be observed that the plasma is found to be much hotter at LHC 
than at RHIC. To study equilibration we have fitted our results with T(r) oc t~", at different 
times. The fitted values of a are given in Table-II and Table-Ill for RHIC and LHC respec- 
tively. Again we see that as time progresses the scaling exponent a increases. As already 
mentioned above, T should not be interpreted as a temperature in the thermodynamic sense 
for T < Teq. This is because the energy- momentum tensor at those early times deviates from 



the ideal fluid form, which is why Eqs. (28-30) do not hold. Nevertheless, T(r) can be taken 
as an indication for the time evolution of the average energy per particle even at r < T^q. 

The evolutions of the number densities are shown in Fig. 4. Not surprisingly, the plasma 
is found to be much denser at LHC than at RHIC. Again, we have fitted our results to 
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n(r) oc r~", at various times. The fitted values of a can be found in Table-II and Table-Ill 
for RHIC and LHC respectively. As already mentioned above, the density nir) decreases 
less fast than oc 1/r. 

The scaling exponents given in Table-II and Table-Ill represent the fitted values of our 
results (e(r), T(r), n(r)) with the functional form r^", at different times. For an equilibrated 
101 dimensionally expanding plasma the scaling exponents are 4/3 for e(r), 1/3 for T(r) 
and 1 for n(r) respectively. It can be seen from Table-II and Table-Ill that our scaling 
exponents approach those values at later times. The reader may judge himself when the 
system is close to equilibrium. Our choice is ~ 4 — 5 fm. 

The time evolutions of the transport cross sections at RHIC and LHC are shown in Fig. 5. 
In our calculation the minijet scale evolves with collision energy, (1.13 and 2.13 GeV at 
RHIC and LHC, respectively; see section ||). Therefore, more energetic partons are produced 
at LHC and the transport cross section is much smaller. These transport cross sections play 
a crucial role in the equilibration of the plasma, cf. Eq. (|TB]). Since the screening mass and 
the momentum scale in the running coupling constant decrease with time the transport cross 
section, in turn, increases strongly. Therefore, despite the decreasing number density the 
relaxation time decreases at later time. In Fig. 6 we have displayed the time evolution of 
the relaxation time rc(r). One observes that the relaxation times at RHIC and LHC do not 
differ by much, despite the much higher density of partons obtained at LHC. 

Finally, in Fig. 7 we present the time evolution of the coupling constant. One observes 
that the coupling seems to be weak at LHC but becomes stronger at RHIC, see also ref. 



One may have to include also higher order processes in the study of equilibration of the 
minijet plasma at RHIC, which will yield faster equilibration as compared to the present 
results. 



VII. SUMMARY AND CONCLUSIONS 

We have studied the production and equilibration of the gluon minijet plasma produced 
in the central region of high-energy nuclear collisions by solving the relativistic Boltzmann 
transport equation. The initial conditions are obtained from pQCD. We have solved the 
transport equation employing a collision term based on 2 — >^ 2 elastic collisions. The coUinear 
divergence in the perturbative cross section to lowest order is removed by incorporating the 
screening of very soft interactions by the medium. The screening mass is calculated from 
the non-equilibrium distribution function of the partons. 

Parton equilibration in high-energy collisions has previously been discussed within the 
Boltzmann equation in refs. [|I^-|2T|. However, ref. made an Ansatz for the time- 



dependence of the relaxation rate, Tc{t) = 9q{t/to)^. It has been found that the results 
depend rather sensitively on the exponent f3. In ||2^ an attempt was made to obtain an 
analytical expression for Tc(r). However, the authors actually use equilibrium distribution 
functions (apart from other approximations) to derive the relaxation time, which is written 
function of temperature. 
In the present work we attempted to reduce such uncertainties by coupling Tc{t) to 
the evolution of the medium itself, leading to self-consistent dynamics. Moreover, Tc{t) is 
determined from the medium by using the number density and transport cross section which 
are computed by using the actual non-equilibrium distribution functions. The treatment of 
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a non-equilibrium medium determining the relaxation time in an explicit and numerically 
computable way is one of the main contributions of our paper. 

Another major difference to previous treatments lies in the choice of the initial time 
and initial conditions. In particular, we have taken into account that according to recent 
arguments regarding parton saturation the minijet scale increases with collision energy, 



which has a significant effect on the initial conditions. In [[19|--|2T[| , for example, the initial 



conditions were taken from the HIJING model at rather late time where it is assumed 
that momentum isotropy is reached (in HIJING the minijet scale Pq is assumed to be energy 
independent). The initial conditions used in the above studies are obtained at Tq= 0.7 fm at 
RHIC and 0.5 fm at LHC. The energy densities at these times are 3.2 GeV/fm^ at RHIC and 
40 GeV/fm^ at LHC. We solve the transport equation starting at tq = 1/po < Tq fm, and do 
not find isotropic momentum distributions at Tq. As our initial conditions at tq = 1/po are 
rather different from those employed in [p!9l-|2l||, our findings regarding equilibration time 
and other bulk properties of the plasma differ considerably from those of the above studies. 

In our study all the quantities such as 'fntij)) ■^(r), ri(r), e(r) are obtained by 

using the non-equilibrium distribution functions. These time dependent quantities are used 
in a self-consistent manner to determine the evolution of the plasma. We have attempted 
to study equilibration without approximations which are valid only in equilibrium, and 
have tried to define most of the quantities in non-equilibrium. For example, we do not 
include a magnetic screening mass in our study because we do not find a closed expression 
for the magnetic screening mass in terms of the non-equilibrium distribution function in the 
literature. Most of the analysis of magnetic screening mass is done in thermal QCD and such 
calculations can not be applied to the initial non-equilibrium stage of high-energy collisions. 
From this point of view we are presently unable to account for magnetic screening in our 
non-equilibrium study. Also, the extensive numerical computations (we solve for the actual 
phase-space distribution function, not only for it's first moment) were presently restricted 
to a momentum independent relaxation time; that is, we considered the static (infinite 
wavelength) limit of the polarization tensor only. 

In the present discussion of thermal equilibration we have restricted ourselves to gg gg 
secondary elastic collisions. We hope to include gg —>■ ggg (and vice versa) inelastic collisions 
in the future, which is necessary to address chemical equilibration. Our present results are 
based on gg gg perturbative elastic collisons of the produced minijets at RHIC and LHC. 

The present study confirms that the plasma is much denser and hotter at LHC en- 
ergy. For central collisions of heavy ions at RHIC and LHC energy the minijet plasma 
appears close to equilibrium at a time Tgg ~ 4 — 5 fm, with a temperature of 220 MeV and 
380 MeV respectively. We do not obtain significantly different kinetic equilibration times 
of the gluon-minijet plasma at RHIC and LHC because the product of comoving parton 
density and transport cross- sect ion, i.e. the kinetic equilibration rate, is similar. However, 
there are differences in other physical quantities such as number density, energy density and 
temperature. Our results indicate somewhat larger equilibration times than those obtained 
in the parton cascade model (with fixed minijet scale po and medium independent cutoff 
for rescattering) and HIJING |2^ . 

For simplicity, we have not incorporated any coherent field in the present study. Including 
a coherent field in the initial condition may decrease the equilibration times and increase 
the energy density, the number density and the temperature of the plasma. We attempt to 
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study the production and equilibration of a QGP at RHIC and LHC with both coherent 
field and incoherent partons taken into account in a forthcoming paper. 
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Figure captions 



FIG. 1. Time evolution of the energy density of the gluon-minijets at RHIC. The sohd 
hne is obtained from the self-consistent solution of the transport equation. The dashed 
line correspond to the free streaming and hydrodynamic limits with the same initial minijet 
conditions at tq = 1/po- 

FIG. 2. Time evolution of the energy density of the gluon-minijcts at LHC. The solid 
line is obtained from the self-consistent solution of the transport equation. The dashed 
line correspond to the free streaming and hydrodynamic limits with the same initial minijet 
conditions at Tq = 1/po- 

FIG. 3. Time evolution of the temperatures of the gluon-minijets at RHIC and LHC. 

FIG. 4. Time evolution of the number densities of the gluon-minijets at RHIC and LHC. 

FIG. 5. Time evolution of the in-medium transport cross sections for the gg — > gg process 
at RHIC and LHC. 

FIG. 6. Time evolution of the kinetic relaxation time Tc(t) at RHIC and LHC. 
FIG. 7. Time evolution of the strong couphng constant ^^(t) at RHIC and LHC. 
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^ (GeV) 


po (GeV) 


A-B no (fm-3) eo (GeV/fm^) 


To (GcV) 


RHIC 200 


1.13 


Au-Au 34.3 


56 


0.535 


LHC 5500 


2.13 


Pb-Pb 321.6 


1110 


1.13 


TABLE L Initial conditions for the pre-equilibrium evolution of the plasma, as 
the number and energy of gluon-minijets with pt > pQ. 


obtained from 


RHIC r (fm) 




a (for e(T)) a (for T(r)) 




a (for n(r)) 


0.5 




1.13 0.286 




0.80 


1.0 




1.15 0.290 




0.76 


1.5 




1.17 0.295 




0.75 


2.0 




1.2 0.302 




0.76 


2.5 




1.22 0.305 




0.77 


3.0 




1.23 0.31 




0.782 


3.5 




1.25 0.312 




0.8 


4.0 




1.26 0.315 




0.83 


4.5 




1.27 0.32 




0.85 


TABLE II. Exponents a 


obtained by fitting our results with the functional form r " at RHIC. 


LHC r (fm) 




a (for e(r)) a (for T(r)) 




a (for n(r)) 


0.5 




1.118 0.29 




0.774 


1.0 




1.164 0.297 




0.75 


1.5 




1.185 0.302 




0.751 


2.0 




1.21 0.306 




0.77 


2.5 




1.227 0.308 




0.8 


3.0 




1.240 0.31 




0.825 


3.5 




1.251 0.314 




0.85 


4.0 




1.262 0.317 




0.875 


4.5 




1.271 0.318 




0.897 


5.0 




1.28 0.32 




0.92 


5.5 




1.28 0.32 




0.92 


6.0 




1.28 0.32 




0.93 


TABLE HI. Exponents a obtained by fitting our results with the functional form 


r-" at LHC. 
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